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^ ABSTRACT 

o 

^ We present a physical characterization of comet 176P/LINEAR, the third member of 

the new class of main-belt comets, which exhibit cometary activity but are dynamically 
^ indistinguishable from main-belt asteroids, to be discovered. Observations show the 

^ object exhibiting a fan-shaped tail for at least one month in late 2005, but then becoming 

inactive in early 2006. During this active period, we measure broadband colors of 
0\ B-V = 0.63 ± 0.02, V - R = 0.35 ± 0.02, and i? - / = 0.31 ± 0.04. Using data from 

when the object was observed to be inactive, we derive best-fit lAU phase function 
Q parameters of H = 15.10 ± 0.05 mag and G = 0.15 ± 0.10, and best-fit linear phase 

^ function parameters of m(l, 1, 0) = 15.35 it 0.05 mag and P = 0.038 it 0.005 mag deg~^. 

Prom this baseline phase function, we find that 176P exhibits a mean photometric 
. ^ excess of ~30% during its active period, implying an approximate total coma dust mass 

of Mrf (7.2 lb 3.6) X 10^ kg. From inactive data obtained in early 2007, we find a 
d rotation period of Prot = 22.23 ± 0.01 hr and a peak-to-trough photometric range of 

Am ~ 0.7 mag. Phasing our photometric data from 176P's 2005 active period to this 
rotation period, we find that the nucleus exhibits a significantly smaller photometric 
range than in 2007 that cannot be accounted for by coma damping effects, and as such, 
are attributed by us to viewing geometry effects. A detailed analysis of these geometric 
effects showed that 176P is likely to be a highly elongated object with an axis ratio 
of 1.8 < b/a < 2.1, an orbital obliquity of e ~ 60°, and a solstice position at a true 
anomaly of i^o = 20 ± 20°. Numerical modeling of 176P's dust emission found that 
its activity can only be reproduced by asymmetric dust emission, such as a cometary 
jet. We find plausible fits to our observations using models assuming ~10 /xm dust 
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particles continuously emitted over the period during which 176P was observed to be 
active, and a jet direction of 180° < aj^t ^ 120° and 5jet ~ —60°. We do not find 
good fits to our observations using models of impulsive dust emission, i.e., what would 
be expected if 176P's activity was an ejecta cloud resulting from an impact into non- 
volatile asteroid regolith. Since for a rotating body, the time-averaged direction of a 
non-equatorial jet is equivalent to the direction of the nearest rotation pole, we find 
an equivalent orbital obliquity of 50° < £ < 75°, consistent with the results of our 
lightcurve analysis. The results of our lightcurve analysis and dust modeling analysis 
are furthermore both consistent with the seasonal heating hypothesis used to explain 
the modulation of 176P's activity. Additional observations are highly encouraged to 
further characterize 176P's active behavior as the object approaches perihelion on 2011 
July 01. 

Subject headings: comets: general — comets: individual (176P/LINEAR = 118401 
(1999 RE70)) — minor planets, asteroids 



1. INTRODUCTION 

The cometary nature of 176P/LINEAR (also known as asteroid 118401; hereafter 176P) was 
discovered on 2005 November 26 as part of the Hawaii Trails Project (HTP; Hsieh et al. 2006; 
Hsieh 2009), a targeted deep-imaging survey aimed at identifying cometary activity in main-belt 
asteroids. Together with two other previously discovered comets — 133P/Elst-Pizarro (hereafter, 
133P) and 238P/Read (hereafter, 238P) — occupying orbits indistinguishable from those of main- 
belt asteroids, the discovery of 176P's cometary nature led to the designation of a new cometary 
class known as main-belt comets (MBCs; Hsieh Sz Jewitt 2006). 

The cometary activity of each MBC is a strong indication of the presence of extant ice (possibly 
preserved in subsurface layers) that has recently become exposed, perhaps by collisions, and is now 
sublimating. Each MBC, however, has a Tisserand parameter (with respect to Jupiter) of Tj > 3, 
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while classical comets have Tj < 3 (Vaghi 1973; Kresak 1980), suggesting that they are unlikely to 
have been captured from the outer solar system and are most likely native to the main belt (Hsieh & 
Jewitt 2006; Jewitt et al. 2009). The limited observational data used to discover the known MBCs 
strongly implies that many more should exist (>100; Hsieh 2009), indicating that present-day ice 
could be widespread in the main asteroid belt. Recent detections of apparent water ice absorption 
at 3.1 /xm in spectroscopic observations of the surface of (24) Themis, the largest asteroid and 
namesake of the Themis asteroid family, to which MBCs 133P and 176P also belong (Hsieh 2009), 
appear to support this conclusion (Rivkin & Emery 2010; Campins et al. 2010), though some believe 
the detected spectral feature could instead be due to non-volatile materials (e.(/.. Beck et al. 2011). 
We note that this uncertainty does not explicitly undermine the case for sublimating ice as the 
driver of MBC activity since, as stated above, that ice is believed to largely reside in subsurface 
reservoirs, and so we would not necessarily expect it to be detectable via reflectance spectroscopy 
of an apparently inactive object like (24) Themis. 



2. OBSERVATIONS 

The 2005 discovery of 176P's comctary nature was made using the 8 m Gemini North telescope 
on Mauna Kea in Hawaii. Following this discovery, confirmation and characterization observations 
were made using Gemini North as well as the University of Hawaii (UH) 2.2 m telescope, also 
on Mauna Kea. Since then, numerous monitoring observations have been made using the UH 
2.2 m and 10 m Keck I telescopes on Mauna Kea, the 3.58 m New Technology Telescope (NTT) 
at the European Southern Observatory (ESO) at La Silla in Chile, and the 2.5 m Isaac Newton 
Telescope (INT) and 4.2 m William Hcrschcl Telescope (WHT) on La Palma in the Canary Islands. 
The object was also observed as part of the HTP using the Lulin 1.0 m telescope in Taiwan just 
one month before its activity was discovered by Gemini, though no clear evidence of activity was 
seen in these observations. All reported observations (detailed further in Table 1) were made in 
photometric conditions. 

Observations with Gemini were made using the imaging mode of the Gemini Multi-Object 
Spectrograph (CMOS; image scale of 0'.'146 pixel~^; Hook et al. 2004), which uses Sloan Digital 
Sky Survey g'r'i'z' filters. UH 2.2 m observations were made using either a Tektronix 2048x2048 
pixel CCD with an image scale of 0'.'219 pixel"^, or the Orthogonal Parallel Transfer Imaging Cam- 
era (OPTIC; 0'.'14 pixcl^^; Tonry et al. 2004). Keck observations were made using the imaging mode 
of the Low- Resolution Imaging Spectrograph (LRIS; Oke et al. 1995), which employs a Tektronix 
2048x2048 CCD (0'.'210 pixer^). NTT observations were made using the ESO Faint Object Spec- 
trograph and Camera (EF0SC2; Buzzoni et al. 1984) which employs a 2048x2048 pixel Loral/Lesser 
CCD (0^'24 pixel"^ using 2x2 binning) behind Bessel BVR broadband filters. INT observations 
were made using the Wide Field Camera which consists of four thinned EEV 2048x4096 pixel 
CCDs (0'.'333 pixer^). WHT observations were made using the Prime Focus Imaging Platform 
(PFIP) which consists of two EEV 2048x4096 pixel CCDs (0'.'24 pixel"^). Lulin observations were 
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made using a VersArray:1300B CCD {Qf!516 pixel"^; Kinoshita et al. 2005). Except where otherwise 
specified, observations were made using standard Kron-Cousins BVRI broadband filters. 

Standard bias subtraction and flat-field reduction was performed on all images. Flat fields were 

constructed from dithered images of the twilight sky, except in the case of our Keck data which used 
flat fields constructed from images of the illuminated interior of the telescope dome. Photometry 
of Landolt (1992) standard stars and field stars was obtained by measuring net fluxes (over sky 
background) within circular apertures, with background sampled from surrounding circular annuli. 
Object photometry was performed using circular apertures of different radii (ranging from 2'.'0 to 
5'.'0), but to avoid the contaminating effects of the coma, background sky statistics were measured 
manually in regions of blank sky near, but not adjacent, to the object. Several (5-10) field stars 
in the comet images were also measured and used to correct for minor extinction variation during 
each night. 



3. RESULTS 

3.1. Active Cometary Behavior 

For all nights on which 176P was observed, we combine individual R- or r'-band images (aligned 
on the object's photocentcr) into deep composite images to assess the level of cometary activity 
present (Fig. 1). Activity is unambiguously detected in observations made of the object between 
UT 2005 November 26 and 2005 December 29. In these observations, no coma is apparent, though 
a fan-shaped tail is clearly visible. No activity in the form of a coma or tail is seen in the next 
available set of observations obtained on 2006 February 03 nor in any monitoring observations made 
since that time. 

Poor seeing during a precovery observation of 176P on 2005 October 24, just one month prior 
to the first confirmed observations of activity, prevents conclusive determination of the object's 

active nature at the time. Photometric evidence (below), however, now suggests that the object 
may have in fact been active. Overall, 176P's behavior is consistent with that of an object that 
exhibits activity near perihelion but not at other points in its orbit (Fig. 2). In all observations 
when 176P is observed to be active, the tail is directed approximately due east, roughly 20° south 
of the anti-solar direction. 

We also measure BVRI colors of 176P on several occasions during its 2005 active period 
(Table 2), finding mean values ofB-V = 0.63 ±0.02, V-R = 0.35 ±0.02, and i?-/ = 0.31 ±0.04, 
and no evidence of substantial variability. This lack of variability is unsurprising because, while no 
significant coma is observed, 176P may very well possess an unresolved coma (c/. 133P; Hsieh et 
al. 2010). If unresolved coma is present, not only will that coma not exhibit color variability as the 
nucleus rotates, it will also dampen any color variability that the underlying nucleus itself might 
exhibit. 
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3.2. Phase Function 

Characterization of an object's photometric dependence on changing solar phase angle, i.e., 
its phase function, permits data obtained at different observing geometries to be compared and 
combined, while also giving insight as to the properties of the surface material of the object itself 
{e.g., Bowell et al. 1989). To compute 176P's phase function, we first consider only photometric 
data obtained when 176P appears to be inactive (i.e., from 2006 Feb 03 onwards; Table 1). 

On several nights, a significant length of time (> 0.6 hr) elapses from the time of our first ob- 
servation of the object to the time of the last observation, during which the object could potentially 
undergo significant rotation. We refer to these data as "lightcurve observations" (Table 3). Obser- 
vations on other nights that do not span a significant amount of time are referred to as "snapshot 
observations" . While on no occasion is 176P seen to reach both the maximum and minimum of its 
rotational lightcurve in a single night of observations, we can nevertheless exploit these lightcurve 
observations to obtain useful information regarding the midpoints of the rotational lightcurve on 
each night. 

For computing 176P's phase function, we aim to use its mean brightness (averaged over a full 
rotation) as a common reference point for data from different epochs. It is obviously not possible 
to know a priori the rotational phase of the object at the time of a given snapshot observation. As 
such, in using photometry from snapshot observations to compute 176P's phase function, we assign 
a maximum uncertainty, a, equal to half of the expected lightcurve range. Am. This allows for 
the range of possible rotational phases up to and including the extreme case where the snapshot 
observation was obtained at either the object's maximum or minimum brightness. 

Unlike in snapshot observations, the rotational phase of lightcurve observations can be con- 
strained by the observed brightness variation, Sm, over the coiurse of the observing period. If 
6m ~ Am, the mean of the minimum and maximum brightness of the object during that set of 
lightcurve observations will be close to the midpoint of the actual lightcurve. We express these 
constraints by assigning uncertainties to the midpoints of our lightcurve observations equal to 
a = (Am — (5m)/2. The final uncertainties listed in Table 3 reflect both these assigned uncertain- 
ties and the uncertainties in the computed mean brightnesses themselves. 

We then perform weighted fits of photometric data obtained from both lightcurve observations 
and snapshot observations to the lAU H,G phase function (Bowell et al. 1989), finding best- 
fit parameters of H = 15.10ib0.05 mag and G = 0.15 it 0.10, which are in agreement with 

the parameters calculated in Hsieh et al. (2009a) with fewer data points and a less sophisticated 
weighting scheme. We also fit our data to a linear phase function, omitting data points obtained 
at solar phase angles of a < 5° at which opposition surge effects are expected, finding best-fit 
parameters of m(l, 1, 0) = 15.35 ± 0.05 mag and (3 = 0.038 ± 0.005, again in agreement with Hsieh 
et al. (2009a). These solutions and the data used to compute them are plotted in Figure 3. 

For all of these calculations, we use Am = 0.70 mag based on the maximum photometric 
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variation observed in lightcurve observations (Table 3). We note, however, that our eventual best- 
fit phase function parameters are not highly sensitive to the value chosen for Am as long as relative 
uncertainties for snapshot and lightcurve observations are preserved. 

3.3. Photometric Confirmation of Activity 

Once a phase function has been computed for 176P using inactive data, we can overlay photo- 
metric data for 176P obtained when the object was active. These data are plotted as solid circles in 
Figure 3 and can be seen to be consistently brighter than the expected mean brightness of the ob- 
ject at those phase angles, confirming the presence of excess flux attributable to cometary activity. 
This type of photometric detection of activity led to the discovery of activity in 95P/Chiron (Bus 
et al. 1988; Tholen et al. 1988; Mcech & Belton 1989; Hartmann et al. 1990) and was also used to 
detect unresolved coma for 133P (Hsieh et al. 2010) and to confirm the reactivation of 238P (Hsieh 
et al. 2011). 

For data taken when 176P was seen to be visibly active (2005 November 26 to 2005 December 
29), we find that the nucleus is, on average, 0.3 =b 0.1 mag brighter than expected, given the H 
and G phase function parameters we derive in §3.2. Following Hsieh et al. (2010), we assume that 
the discrepancy between the observed magnitude, ruohs^ and the expected magnitude, ffiexpi is due 
to dust contamination, and that the albedos of the nucleus and dust are equal. We then calculate 
that the scattering surface area of the dust coma (and any portion of the tail also contained within 
the photometry aperture), A^, is ~ 30% that of the inactive nucleus, An, using 

An 

Using An = vrrg = (1.3 ib 0.2) x 10^ m^, where Ve = 2.0 ± 0.2 km is the effective radius of 176P's 
nucleus (Hsieh et al. 2009a), we therefore find A^ji = (0.4 it 0.2) x lO'^ m^. Assuming a typical dust 
grain radius of 10 //m — similar to dust grain radii found for 133P (Hsieh et al. 2004) and consistent 
with dust modeling results for 176P (§4) — and a bulk grain density of = 1300 kg m~^ (Hsieh et 
al. 2004), we find an approximate total coma dust mass of ~ (7.0 it 3.5) x 10^ kg, comparable 
to coma dust masses found for 133P (c/. Hsieh et al. 2010). 

3.4. Rotation 

3.4-1. Period Determination 

Given the minimal nucleus-obscuring coma observed during 176P's 2005 active phase, we 
performed phase dispersion minimization analysis {e.g., Stellingwerf 1978) on the photometry data 
we obtained during that time in an attempt to ascertain the nucleus's synodic rotation period, but 
were unable to identify a period that would allow us to produce a convincing lightcurve. At the 
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time, we attributed this failure to find a plausible rotation period to a probable extremely slow 

rotation rate, which made it difficult to sample a significant portion of the rotational lightcurve 
in a single night (where the object was visible for ~4 hr in 2005 December), and possible aliasing 
caused by a rotation period close to the Earth's own daily 24-hr rotation period. 

A more suitable data set for determining 176P's rotation period was obtained on two nights 
in 2007 February and two additional nights in 2007 March. In addition to affording longer nightly 
visibility windows (~7 hr in February and ~5 hr in March) , the timing of the observations over two 
consecutive months allowed us to avoid both aliasing effects and complications due to significant 
changes in viewing geometry. Applying the same phase dispersion minimization techniques as before 
and phasing our data to candidate rotation periods to assess plausibility, we find a likely rotation 
period of Prot = 22.23 it 0.01 hr (assuming the lightcurve is double-peaked) and a peak-to-trough 
photometric range of Am ~ 0.7 mag (Fig. 4; where data for 2007 May 19 is overplotted to check for 
consistency) . We estimate the uncertainty on the period by slowly varying the optimum period and 
re-phasing the data until we determine that the lightcurve is noticeably out of phase. We also find 
possible but less likely candidate periods of Prot = 22.57 hr, Prot = 22.83 hr, and Prot = 23.20 hr, 
all of which produce lightcurves with larger phase dispersions than for Prot = 22.23 hr. 

The lightcurve we find using our most likely rotation period is incompletely sampled, but 
otherwise appears convincing. These results confirm our preliminary findings from the analysis 
of our 2005 data that the rotation of the object is extremely slow, but also show that the object 
is relatively elongated, with an axis ratio of a/b > lO^-^^"*, or a/b > 1.9. Interestingly, the 
lightcurve also exhibits features such as a V-shaped minimum and an inverted U-shaped maximum 
that are qualitatively like those of contact binaries {e.g., Sheppard & Jewitt 2004; Mann et al. 
2007; Lacerda 2008), though the photometric range we observe for 176P is not quite as large as is 
normally expected for such systems (Am > 0.9 mag; Lacerda 2008). 

3.4.2. 2005 Photometry Revisited 

Having achieved a reasonable lightcurve using 2007 data (§3.4), we turn back to our 2005 data 
and attempt to phase them using the same rotation period (Fig. 5). When phased together, the 
data fit together far less well than the 2007 data, but the resemblance to the 2007 lightcurve is 
clear. We attribute much of the scatter in the data to the extreme sensitivity of photometry of 
active objects to seeing, where even small changes in image quality from exposure to exposure 
can cause significantly different amounts of dust contamination to be present in each photometry 
aperture. We believe this effect to be the cause of the slight brightness excess toward the second 
half of our 2005 December 25 observations (Fig. 5). On 2005 December 24, we find an overall 
brightness excess for the comet in the range of ~ 0.1 — 0.2 mag (see Table 1), but given that the 
comet was passing through a particularly dense star field on this night, we attribute this excess 
to contamination from unseen faint field stars, and omit this data from Figure 5 for clarity. An 
active cometary object may of course be expected to occasionally exhibit short-lived outbursts, but 
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given the random (non-correlated with rotational phase) and short-lived nature of the observed 
photometric anomalies, we find that they arc adequately explained by fluctuations in the nightly 
seeing and, on occasion, by faint field star contamination. 

Significantly, the phased 2005 data appear to exhibit a far smaller photometric range than the 
2007 data, showing a peak-to-trough variation of Am 0.2 mag, down from Am ~ 0.7 mag in 
2007. As the observed photometric range of a rotating nucleus can be damped by a superimposed 
coma of constant brightness (determined to be present for 176P in §3.3), before evaluating the 
physical implications of this change in observed photometric range, we must first determine how 
much of this effect could be due to coma damping. 

Noting that the measured flux from an active comet. Fobs, consists of the sum of the fluxes 
from the nucleus and dust contained within the photometry aperture, we can write 

-l0.4Amoi,s Fobg ij^ax ( Fn^max 



-|^QU.4/lmoi,s _ - oos,max _ i ^ n,max > ^ a i 
Fohs,min \ Fji^min + FdJ 
where Aniobs is the observed photometric range of the lightcurve, Fobs,max and Fobs,min are the 
maximum and minimum fluxes, respectively, observed for the active nucleus, Fn^max and Fn^min are 
the maximum and minimum fluxes, respectively, for which the nucleus is responsible, and the flux 
due to dust, F^, is assumed to be constant. 

For a true nucleus axis ratio (i.e., uncontaminated by coma) of {a/b)n, nucleus fluxes can be 
expressed as follows: 

F -('^\"\ (^) 
^n,max — \ i I n 



n 



and 

Fn,min = ( T ) i^) 



.b 

where F„ is the average flux of the nucleus over one full rotation. 
Then, substituting Equations 3 and 4 into Equation 2, we flnd 



Q + ^ (1 _ 100-4Am„,.) _ 100.4Am„,. ^ (5) 

1/2 

which can then be solved for {a/h)n using standard techniques for solving quadratic polynomials, 
where F^/Fn = Ad/ An is given by Equation 1. Taking the positive root of Equation 5 (since Equa- 
tions 3 and 4 must produce positive values), we flnd {a/b)n = 1.27, meaning that even accounting 
for coma contamination, the photometric range of the bare nucleus in 2005 is Am„ = 0.26, still 
much smaller than that observed in 2007. 



'a\i/2 



3.4-3. Implications for Pole Orientation and Physical Nature 



The plausibility of the seasonal heating hypothesis as a mechanism for modulating MBC ac- 
tivity (flrst proposed for 133P; Hsieh et al. 2004), depends crucially on each object having a pole 
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orientation compatible with its activity profile. Pole orientations are currently unknown for all 
MBCs, though some constraints placed on 133P's pole orientation by Toth (2006) are thus far 
consistent with seasonal modulation of that MBCs activity. 

In the case of 176P, we can use the rotational lightcurves shown in Figures 4 and 5 to place con- 
straints on its pole orientation. These figures show that the photometric range of 176P's lightcurve 
changed appreciably between 2005 (Am ~ 0.2 mag) and 2007 (Am ~ 0.7 mag), which suggests 
that the object has significant obliquity. The available data (incomplete lightcurves at two epochs) 
are not sufficient to uniquely determine the spin pole direction of 176P. For example, we are unable 
to break the degeneracy between north and south pole. However, we can still use those data to 
obtain useful constraints on the rotational pole of 176P. 

In particular, we are interested in constraining the shape and obliquity of 176P and verifying 
if the data are consistent with the seasonal hypothesis for MBC activity. According to the seasonal 
hypothesis, 176P should have non-zero obliquity and be active close to one of its solstices. Wc thus 
elect to parameterize the orientation of the spin pole of 176P by its obliquity, £, and by the true 
anomaly of the solstice, uq, which we assume to be the northern hemisphere summer solstice. For 
instance, if 176P were at solstice exactly at perihelion, then z/q = 0°. We assume the already low 
orbital inclination {i = 0.23°) of 176P to be exactly zero for the purpose of this analysis. 

To investigate the lightcurve behavior of 176P we employ the simulations described in Lacerda 
&L Jewitt (2007). These simulations place triaxial ellipsoids (semi-axes a > b > c) at pre-selected 
illumination and viewing angles and register the integrated reflected flux as a function of rotational 
phase (assuming simple principal axis rotation about the ellipsoids' short axes) to extract modeled 
lightcurves. In addition to observing geometry (completely defined by e, uq, and the phase angle 
a), the model takes as parameters the axis ratios (0 < c/a < b/a < 1) of the triaxial ellipsoid 
and the scattering law. For simplicity, we use prolate ellipsoids (c/a = b/a) and a Lommel-Seeliger 
"lunar" scattering function which has no free parameters and is appropriate for simulating the low 
albedo {pji = 0.06 ± 0.02; Hsieh et al. 2009a) surface of 176P. In what follows, we compare the 
model lightcurves with the data (Figs. 4 and 5) from 2005 (chiefly taken at point (c) in Fig. 2) and 
2007 (taken at points (h), (i) and (j) in Fig. 2). 

Due to the underconstrained nature of this problem given the available data, we only consider 
two limiting scenarios, one in which we simply assume that the lightcurve amplitude observed 
in 2007 represents the object's maximum range, placing it at equinox during those observations, 
and a second scenario that is physically motivated where we assume that the object was receiving 
maximal heating at the time when our 2005 observations showed it to be active, placing solstice at 
perihelion. 

First, we consider the possibility that 176P was close to equinox during the 2007 measure- 
ments, implying that the solstice occurs around true anomaly = 21°, i.e., just after the 2005 
observations. In that case we can use the photometric range in 2007 to obtain a direct estimate of 
the b/a axis ratio from the relation Am = — 2.5 log(6/a). Using Am ~ 0.7 mag, we find b/a ~ 0.52. 
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Figure 6 (left panel) shows the 2007 lightcurve data overplotted on three simulations spanning 
0.54 < b/a < 0.58. The reason why the best-fit model shapes are slightly less elongated than the 
purely geometric estimate, b/a = 0.52, is that limb darkening effects contribute to increase the 
photometric variability beyond the geometric expectation. Each model was calculated using the 
mean phase angle and true anomaly of the 2007 observations. Under these conditions, the obliquity 
is approximately determined by the 2005 "solstice" data. An obliquity of e = 0° would result in a 
nearly constant photometric range throughout the orbit, whereas e = 90° would imply negligible 
photometric variation in 2005. We find that an obliquity e ~ 60° for an object with an axis ratio 
b/a - 0.56 ± 0.03 fits the 2005 data weU (Fig. 6, right panel). 

We now consider a second scenario in which solstice occurs at perihelion, z/q = 0°, just before 
the period of activity in 2005. We reason that if the activity is powered by the seasonal heating of 
sub-surface ice, then it should begin at or shortly after the time of maximum localized heating, i.e., 
the solstice. If uq = 0°, then equinox takes place at true anomaly u = 90°, or about 20° before the 
2007 data were taken. In this scenario, 176P must be more elongated than determined above. In 
Figure 7 we show this to be the case: to simultaneously fit the 2005 and 2007 data, for an obliquity 
of e = 60°, an axis ratio 6/o 0.52 is required. 

In both scenarios, we note that models using significantly smaller obliquities {e.g., e = 45°) 
cannot simultaneously fit the 2005 and 2007 {e.g., Fig. 8). Higher obliquities are possible, but 
require a more elongated nucleus shape for 176P and a later solstice location along the orbit. 
Figure 9 shows the result of simulations for solstice at fo = 40°, obliquity e = 85°, and shapes 
0.33 < b/a < 0.48. However, objects with such extreme axis ratios (3:1, 2.5:1 and 2.1:1 in Fig. 9) are 
hydrostatically unstable (Jeans 1919) and require significant material strength to retain their shape. 
Alternatively, such extreme configurations may be explained by contact binaries (Lacerda & Jewitt 
2007), but our current data set is insufficient for confirming or ruling out this possibility. Solstice 
positions significantly outside the interval 0° < I'o ^ 40° (or equivalently, 180° < fo < 220°, due to 
the fact that "summer" and "winter" solstices are indistinguishable in this analysis) are untenable 
as the 2005 data with its smaller Am would fall as close or closer to the equinox than the 2007 
data with its larger Am, when in fact Am should reach a maximum at equinox. 

In conclusion, our data strongly suggest that 176P is highly elongated, close to a 2:1 axis ratio, 
and has significant obliquity, close to £ = 60°. The solstice positions that best fit the data are 
located close to the active portion of the orbit, around fo = 20 ± 20°, placing the other solstice at 
viso = 200 ± 20°, and equinoxes at ugo = 110 ± 20° and 1^270 = 290 ± 20°. These constraints on the 
obliquity and solstice position of 176P are consistent with the seasonal hypothesis to explain the 
activity of MBCs. 
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4. 



DUST MODELING 



In order to place quantitative constraints on 176P's dust emission, we generate a series of 
numerical models for the dust and attempt to match them to our observations. Due to the limited 
observational data available and the resulting underconstrained nature of this dust modeling effort, 
we recognize from the outset that it will not be possible to achieve an exact model of 176P's 
activity. As such, we purposefully formulate our modeling strategy to simply achieve constraints 
on particular key properties such as grain sizes, ejection velocities, and the temporal behavior of 
the dust emission. 

Standard dust modeling (e.g., Finson & Probstein 1968) typically makes use of syndyne curves, 
which are lines representing the positions of particles of constant sizes ejected at different times 
where ejection velocities are assumed to be zero, and synchrone curves, which are lines representing 
the positions of particles of different sizes ejected at the same time with zero velocity. In the case of 
176P, an extremely low inclination means that all synchrone and syndyne curves actually overlap 
with one another in the object's orbit plane, making it very difficult to use this type of analysis to 
study 176P's dust emission. We note, however, that the orientation of 176P's dust tail does not 
coincide with its orbit plane, and thus, does not actually coincide with any syndyne or synchrone 
curves (Fig. lb,c,d), and therefore we must consider asymmetric dust emission (as in a jet), similar 
to that considered in our previous analysis of 238P (Hsieh et al. 2009b). The methods used in 
the following analysis are the same as those used in that previous work except where otherwise 
specified. 

Assuming jet-driven dust emission, we assume that the central axis of the cone points in a 
particular direction in the inertial frame (i.e., toward right ascension aj^t, and declination Sjet)- 
For a rotating body and a non-equatorial jet, it should be noted that the effective time-averaged 
jet position in this model actually simply corresponds to the nearest rotational pole, where the 
derived jet width will be somewhat larger than the actual jet width due to the sweeping motion 
caused by the body's rotation (c/. Fig. 10). For simplicity, we assume that Earth-bound observers 
are situated outside the cone defined by the rotating jet, and that the strength of the jet does not 
vary with rotational phase. 

The time-averaged jet is presumed to originate at the subsolar point and have a half-opening 
angle of w = 45°. We furthermore assume that dust particles are released homogeneously from a 
spherical body and that emission only occurs on the day-time side of the body. Dust particle sizes 
are parameterized in standard fashion where /3 denotes the ratio of a particle's acceleration due to 
solar radiation pressure to its acceleration due to gravity. Terminal velocities, Vej, of ejected dust 
particles are given by 



where rh is the heliocentric distance, vq is the reference ejection velocity in m (assumed here 
to be vo = 25 m s~^ based on our work with 238P; Hsieh et al. 2009b) of particles with /5 = 1 
at r/i = 1 AU, and ui and U2 (assumed here to be 1x1 = ^2 = 0.50) are the power indices of the 




(6) 
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reference ejection velocity dependence on ^ and Vh- We use an exponential size distribution with 

an index of g = —3.5, minimum /3 value of /3mm = 9 x 10~^, and 

and an exponential dust production rate dependence on heliocentric distance with an index of 
k = 2>. Dust emission is assumed to begin 1 year prior to perihelion passage and is terminated on 
2006 January 1 (three days after 176P was last observed to be active). We assume such an early 
start date to allow time for the dust to evolve to a "steady" state, but given the short dissipation 
timescalcs for micron-scale dust particles, also note that most of the dust that is actually observed 
is likely to have been ejected far more recently. 

We show resulting modeled dust clouds for 2005 November 26 and 2005 December 29 for 
different jet orientations in Figures 12 and 13, respectively. Comparing these to the observations in 
Figure 11, we see that reasonable matches to our data are provided by models with jet directions 
of 180° < ctjet ^ 120° and 5jet ~ —60°, or approximately oriented towards the Sun during the 
comet's active period. Crucially, all model scenarios also reproduce the non-detection of activity 
on 2006 February 03 (Fig. 1), where all dust particles in each scenario have dispersed from the 
field of view by that time. This rapid disappearance of dust activity did not occur from a previous 
set of test scenarios for which we assumed much larger particles were present {(imin = 10""^). We 
thus conclude that the dust particles ejected by 176P must have been on the order of ^■lO ;um 
in size (corresponding to /3 ~ 10^^), similar to the dominant particle sizes found for 133P (Hsieh 
et al. 2004). Using Equation 6, we find an approximate ejection velocity for these particles of 
Vej ~ 5 m s"-*^, or somewhat faster than the Vgj ~ 1 m s"-*^ ejection velocities found for 133P. 

To determine an approximate dust production rate, we refer back to our estimate for 176P's 
unresolved coma mass in §3.3 of ~ 7 x 10^ kg. This mass estimate is based on nucleus 
photometry employing 3'.'0-radius photometry apertures, which at the distance of 176P in December 
2005, correspond to physical radii of Vcoma ^ 5000 km. The time that it would take a dust particle 
to cross the photometry aperture can then be estimated as tcross = fcoma/vej ~ 1 x 10^ s, or about 
12 days. Using this as the approximate timescale on which dust in the coma must be replaced by 
new dust production from the nucleus, we therefore find an approximate dust production rate of 
cross 0.07 kg s ^ . This production rate is similar to that found for 133P by (Hsieh 
et al. 2004), but is also likewise probably accurate, at best, to an order of magnitude. 

We are further interested in whether 176P's dust emission is only consistent with continuous 
emission, i.e., that would be expected if the activity is driven by the sublimation of ice, or whether 
it can also be reproduced by an impulsive emission event, i.e.. if 176P's dust tail solely consists 
of ejecta produced in an impact by another asteroid. To test this, we produce another series of 
models in which dust emission is limited to a single burst of particles. We test three different size 
distributions — 10"^ < P < 10-^ 10"^ < 13 < 10-\ 10"^ < ^ < 10"^ — where dust in each 
model is ejected on a single day (2005 November 15), which is four weeks after perihelion and 11 
days before activity was first observed (Fig. 14). We then continue to follow the evolution of the 
dust and compare its appearance in each model to observations on 2005 December 29 and 2006 
February 03. 
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We clearly see that for the models including large particles {Pmin = 10~ and /3min = 10~ ; 
second and third columns of Fig. 14), the appearance of the comet remains approximately constant 
between November and December, as is observed. However, the comet retains its appearance 
through February as well (consistent with our continuous emission models), which is inconsistent 
with observations. If only smaller particles are included {Pmin = 10~^; fourth column of Fig. 14), 
we note that the dust cloud does dissipate appreciably by February, but also undergoes significant 
dissipation by December, inconsistent with observations. Thus, we find that a single burst of 
particles cannot produce a dust cloud that simultaneously remains constant between 2005 November 
26 and 2005 December 29, but then dissipates by 2006 February 03. 

We therefore conclude that 176P's activity is most likely driven by an extended emission event, 
in which 176P's dust tail is continuously replenished by small, fast-dissipating particles between 
2005 November and December, allowing it to maintain its appearance, but where dust emission 
ceases sometime in December, leading to the nearly complete dissipation of the dust cloud by the 
time the comet is observed again on 2006 February 03. This behavior strongly suggests that 176P's 
activity is driven by the sublimation of volatile ices, as was previously found for fellow MBCs, 133P 
and 238P (Hsieh et al. 2004, 2009b). If this conclusion is correct, we would also expect 176P to 
exhibit repeated activity, similar to 133P and 238P (Hsieh et al. 2010; Hsieh 2010), perhaps during 
its upcoming perihelion passage (2011 July 01). Observations of renewed activity at this time will 
provide strong support for sublimation at the cause of 176P's activity. Conversely, however, if 
no activity is observed, reconsideration of impact-generated dust emission as the source of 176P's 
activity may be necessary. 

In terms of pole orientation, our dust modeling analysis is remarkably consistent with the 
results of our lightcurve analysis (§3.4.3), as a jet orientation of 180° < ajet ^ 120° and ^jet ~ —60° 
that is effectively coincident with the object's rotational pole corresponds to an orbital obliquity 
of 50° ^ e ^ 75°. The consistency of these results is encouraging support for their accuracy, and 
strongly suggests that 176P is in fact a high-obliquity object that is at solstice around the time that 
a localized site in the "summer-time" hemisphere experiences active sublimation, just as described 
under the seasonal heating hypothesis for MBCs. 

5. DISCUSSION 

5.1. Activity Profile 

Little is currently known about the activity profile of 176P. During 176P's 2005 active period, 
dust emission appeared to end between v ~ 19° and v ~ 28° in true anomaly. No conclusive 
evidence of activity was observed from Lulin Observatory when the comet was at a true anomaly 
of = 1.4°, though a photometric analysis (§3.3) suggests that activity may have been present but 
simply escaped detection, perhaps due to poor seeing. As such, the point at which 176P's activity 
began is essentially unconstrained. If 176P is similar to 133P and exhibits recurrent, seasonally- 
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modulated dust emission that persists over ~ 90° of its orbit, though, we might have expected it to 

resume cometary activity as early as 2010 September when it reached v ~ 290°. If, however, 176P 
is active over a much smaller portion of its orbit than 133P (as it might, for example, if shadowing 
due to local topography also has a strong activity-modulating effect), then activity may not return 
until closer to the object's perihelion passage on 2011 July 01. Continued monitoring of 176P is 
highly encouraged to help clarify this issue. 

5.2. Comparison to other MBCs 

The activity of 176P is similar to 133P (Hsieh et al. 2004, 2010) in terms of dust output (c./. 

§4), but much weaker than highly active MBCs like 238P (Hsich et al. 2009b, 2011), P/2008 Rl 
(Garradd) (hereafter P/Garradd; Jewitt et al. 2009), or P/2010 R2 (La Sagra) (hereafter P/La 
Sagra). Recently observed comet-like activity for two other objects in the main asteroid belt — 
P/2010 A2 (LINEAR) and (596) Scheila (Birtwhistle et al. 2010; Larson 2010) — is likely to be due 
to impact-generated ejecta clouds, and not sublimation-driven dust emission (Jewitt et al. 2010; 
Snodgrass et al. 2010; Jehin et al. 2011), and so we omit these so-called "disrupted asteroids" from 
this discussion. 

Differences in MBC activity strength could be due to low-activity MBCs actually being in- 
trinsically less icy than highly active MBCs, or perhaps simply activated less recently than highly 
active MBCs, resulting in greater depletion of their current active sites. Observations of highly 
active MBCs during successive future active periods to track their decline in strength, if any, will 
help to clarify this issue. Noting the distribution of activity levels of MBCs discovered in the future 
may also be informative. A high ratio of low-activity objects to highly active objects could indicate 
asymptotic declines in activity strength as individual active sites become depleted. On the other 
hand, a flatter distribution of MBCs of varying strengths, or evidence that activity strength has 
a particular spatial dependence, could indicate that activity strength is more strongly dependent 
on the intrinsic ice content of each body (and by extension, the ice content of the surrounding 
population) . Any analysis of this nature will also have to account for observational discovery biases 
towards more highly active objects. 

The activity profile discussed above (§5.1) is particularly interesting in relation to 133P in that 
while 133P becomes active shortly before perihelion, 176P becomes inactive shortly after perihelion 
(Fig. 2). This fact is significant because if 176P's thermal inertia is low, it should reach its maximum 
overall temperature close to perihelion. If thermal inertia is high, 176P's temperature may not be 
at its maximum at perihelion, but should instead be increasing at this point of close approach to 
the Sun. In either case, the fact that 176P's activity ends (or at least declines below detectable 
levels) shortly after the object's overall temperature has peaked or is even still increasing strongly 
favors a scenario in which activity is modulated by seasonal effects, as proposed for 133P (Hsieh et 
al. 2004, 2010), rather than proximity to the Sun and overall temperature. 
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One caveat, though, is that we do not actually know whether the duration of 176P's activity 
is similar to 133P's (i.e., persisting over ~25% of its orbit). The active portion of 176P's orbit 
might still be centered on perihelion if, in addition to ending just after perihelion, it also only 
starts just before perihelion, giving it a much shorter active arc than 133P. Such a shorter active 
arc could perhaps be due to a much smaller (or older and therefore more depleted) supply of 
exposed volatile material, or topography imparting much more local shadowing, and therefore a 
much sharper seasonal effect, than on 133P. Again, as discussed above, continued monitoring of the 
object in search of resumed activity as it approaches its next perihelion passage will be essential 
for resolving these questions. 

6. SUMMARY 

We present a physical analysis of main-belt comet 176P/LINEAR based on optical observations 
and numerical modeling, finding the following key results: 

1. In optical imaging data obtained between late 2005 and mid-2009, we detect the presence of 
cometary activity for 176P between 2005 November 26 and 2005 December 29, and do not 
detect activity between 2006 February 03 and 2009 May 03, placing the apparent cessation 
point for 176P's activity at 20°-25° past perihelion. The turn-on point for activity remains 
unknown at this time. 

2. Using photometric data obtained when no activity was detected for 176P, we find best- 
fit lAU phase function parameter values oi H = 15.10 it 0.05 mag and G = 0.15 it 0.10, 
and best-fit linear phase function parameter values of m(l,l,0) = 15.35 =b 0.05 mag and 
13 = 0.038 ± 0.005 mag deg-^ 

3. Using data obtained in 2007 when no activity was detected, we also find a rotation period solu- 
tion of Prot = 22.23±0.01 hr and a photometric range of Am ~ 0.7 mag. Using this rotational 
period to phase data from 2005 when 176P was active, we find a much smaller photometric 
range of Am ~ 0.2 mag which cannot be fully accounted for by coma contamination and as 
such is attributed by us to viewing geometry effects. From these lightcurve measurements, 
we derive a likely orbital obliquity of £ '-^ 60° and an axis ratio of 1.8 < h/a < 2.1. 

4. In performing dust modeling of 176P's activity, we noted that the orientation of the dust tail 
does not coincide with any syndyne or synchrone curves, requiring us to consider asymmetric 
dust emission, e.g., a cometary jet. We find plausible fits to our data using a continuous dust 
emission model with a particle size range of 9 x 10~^ < j3 < 10~^ and a reference ejection 
velocity of = 60 m s~^. Our data are not consistent with any impulsive dust emission 
scenarios that we investigated. We find a likely jet orientation (and therefore rotational pole 
orientation) of 180° < ajet 120°and 5jet ~ —60°, in good agreement with our lightcurve 
analysis above. 
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Table 1. Observation Log 



UT Date Tel.'' t" Filters Active? u" As Opj' mflJ mji(l,l,a)^ 



2005 Oct 18 


Perihelion 












0.0 


2.58 


1.59 


1.5 


-0.1 






2005 Oct 24 


Lulin 


10 


3000 


R 


1.5 


yes? 


1.4 


2.58 


1.60 


4.3 


-0.1 


18.20±0.01 


15.12±0.01 


2005 Nov 26 


Gemini 


2 


240 


r' 


0.6 


yes 


10.1 


2.59 


1.82 


16.4 


-0.1 


19.11±0.04 


15.74±0.04 


2005 Dec 22 


UH/Tek 


7 


2100 


R 


1.1 


yes 


16.8 


2.60 


2.12 


21.1 


-0.1 


19.65±0.01 


15.95±0.01 


2005 Dec 24 


UH/Tek 


26 


7800 


R 


1.0 


yes 


17.3 


2.60 


2.15 


21.3 


-0.1 


19.46±0.01 


15.73±0.01 


2005 Dec 25 


un^ ±eK 




yyuu 






yes 


1 7 fi 


9 fin 


2 16 


91 A 


— n 1 

— U.J. 


iQ ao-i-ci ni 


J-O.ooinu.u J. 


2005 Dec 26 


UH/Tek 


31 


9300 


BVRI 


1.0 


yes 


17.8 


2.60 


2.17 


21.5 


-0.1 


19.62±0.01 


15.86±0.01 


2005 Dec 27 


UH/Tck 


29 


8700 


R 


0.8 


yes 


18.1 


2.60 


2.19 


21.5 


-0.1 


19.59±0.01 


15.81±0.01 


2005 Dec 29 


Gemini 


34 


3060 


a'r'i'z' 


0.7 


yes 


18.6 


2.60 


2.21 


21.7 


-0.1 


19 62+0 01 


1 5 SO+0 01 


2006 Feb 03 


UH/Tck 


11 


3300 


R 


1.1 


no 


27.7 


2.63 


2.71 


21.2 


0.0 


20.25±0.01 


15.98±0.01 


2006 Aug 31 


UH/Tek 


2 


600 


R 


0.8 


no 


75.1 


2.93 


3.55 


14.3 


0.0 


21.23±0.08 


16.14±0.08 


2006 Sep 02 


UH/Tek 


2 


600 


R 


0.9 


no 


75.5 


2.94 


3.53 


14.6 


0.0 


21.09±0.05 


16.00±0.05 


2006 Dec 11 


UH/Opt 


30 


9000 


R 


1.0 


no 


94.4 


3.12 


2.42 


14.5 


-0.1 


20.01±0.01 


15.62±0.01 


2006 Dec 16 


UH/Opt 


11 


3300 


R 


0.9 


no 


95.3 


3.13 


2.38 


13.3 


-0.1 


20.13±0.01 


15.77±0.01 


2006 Dec 18 


UH/Opt 


3 


900 


R 


1.3 


no 


95.6 


3.14 


2.36 


12.8 


-0.1 


20.09±0.09 


15.74±0.09 


2007 Jan 27 


Keck 


3 


240 


R 


0.9 


no 


102.5 


3.21 


2.23 


0.8 


-0.1 


19.50±0.01 


15.23±0.01 


2007 Feb 15 


UH/Tek 


30 


9000 


R 


0.9 


no 


105.6 


3.25 


2.33 


7.5 


0.0 


19.87±0.01 


15.47±0.01 


2007 Feb 16 


UH/Tek 


49 


14700 


R 


1.2 


no 


105.8 


3.25 


2.33 


7.8 


0.0 


19.93±0.01 


15.53±0.01 


2007 Mar 21 


UH/Tek 


35 


10500 


R 


0.8 


no 


111.1 


3.31 


2.72 


15.4 


0.0 


20.71±0.01 


15.94±0.01 


2007 Mar 22 


UH/Tck 


51 


15300 


R 


1.3 


no 


111.3 


3.31 


2.73 


15.5 


0.0 


20.82±0.01 


16.04±0.01 


2007 May 19 


UH/Tek 


7 


2100 


R 


1.1 


no 


120.2 


3.41 


3.64 


16.1 


0.1 


21.57±0.05 


16.10±0.05 


2008 Jun 29 


NTT 


2 


360 


R 


1.3 


no 


173.2 


3.80 


3.80 


15.4 


0.1 


21.68±0.07 


15.88±0.07 


2008 Jun 30 


NTT 


3 


540 


R 


1.0 


no 


173.3 


3.80 


3.81 


15.3 


0.1 


21.70±0.05 


15.90±0.05 


2008 Jul 01 


NTT 


3 


540 


R 


1.1 


no 


173.4 


3.80 


3.83 


15.3 


0.1 


21.63±0.05 


15.81±0.05 


2008 Aug 25 

2009 Jan 23 


Aphelion .. 












180.0 


3.81 


4.52 


9.9 


0.0 






WHT 


4 


240 


R 


0.8 


no 


198.0 


3.77 


4.01 


14.1 


0.0 


21.47±0.10 


15.57±0.10 


2009 May 03 


INT 


2 


600 


R 


1.8 


no 


210.5 


3.69 


2.70 


3.8 


0.1 


20.30±0.04 


15.31±0.04 


2011 Jul 01 


Perihelion 












0.0 


2.58 


2.98 


19.4 


0.1 







^Telescope used 

''Number of images 

'^Total effective exposure time 

■^FWHM seeing in arcsec 

''True anomaly in degrees 

^Heliocentric distance in AU 

^Geocentric distance in AU 

''Solaj: phase angle (Sun-176P-Earth) in degrees 

'Orbit plane angle (between the observer and object orbit plane as seen from the object) in degrees 
JMcan (in magnitude space) of maximum and minimum ij-band magnitudes measured for nucleus 

*^Inferred reduced i?-band magnitude (normalized to iJ = A = 1 AU) at midpoint of full photometric range (assumed to be 0.70 mag) of 
rotational lightcurve 
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Table 2. Color Measurements 



UTDate B-V V-R R-I 



2005 Dec 25 
2005 Dec 25 
2005 Dec 26 
2005 Dec 29 
Mean 



0.61±0.03 
0.61±0.03 
0.64±0.03 
0.64±0.02 



0.32±0.02 
0.36±0.03 
0.37±0.02 
0.36±0.02 



0.27±0.03 
0.27±0.04 
0.35±0.03 
0.32±0.02 



0.63±0.02 0.35±0.02 0.31±0.04 



Table 3. Lightcurve Observations 



l^T Dale 










2006 Feb 03 


1.65 


0.13±0.06 


20.28±0.06 


16.02±0.29 


2006 Dec 11 


4.23 


0.39±0.12 


20.02±0.12 


15.63±0.20 


2006 Dec 16 


1.98 


0.28±0.05 


20.10±0.05 


15.74±0.22 


2007 Feb 15 


6.71 


0.51±0.04 


19.92±0.04 


15.62±0.04 


2007 Feb 16 


7.34 


0.55±0.06 


20.01±0.06 


15.68±0.06 


2007 Mar 21 


4.66 


0.63±0.07 


20.77±0.07 


16.00±0.08 


2007 Mar 22 


5.12 


0.55±0.10 


20.83±0.10 


16.05±0.13 


2007 May 19 


0.86 


0.31±0.22 


21.53±0.22 


16.06±0.32 



^Time spanned by observations (hr) 

^Photometric range between maximum and minimum i?-band 
magnitudes measured for the nucleus 

'^Midpoint between maximum and minimum i?-band magni- 
tudes measured for the nucleus 

^Inferred reduced i?-band magnitude (normalized to i? = A = 
1 AU) at midpoint of full photometric range (assumed to be 0.70 
mag) of rotational lightcurve 
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Fig. 1. — Composite images of 176P from i?-band images taken during observations detailed in Table 1. 
Each image is 0^5 x 0^5 with 176P at the centre, with arrows indicating north (N), east (E), the negative 
heliocentric velocity vector {—v), and the direction towards the Sun (0). Images shown comprise (a) 3000 s 
of exposure time on the Lulin 1.0 m telescope, (b) 240 s on the 8 m Gemini North telescope, (c) 37500 s 
on the University of Hawaii 2.2 m telescope, (d) 2520 s on Gemini North, (e) 3600 s on the UH 2.2 m, (f) 
2400 s on the UH 2.2 m, (g) 9000 s on the UH 2.2 m, (h) 240 s on the 10 m Keck I telescope, (i) 24600 s 
on the UH 2.2 m, (j) 25800 s on the UH 2.2 m, (k) 3000 s on the UH 2.2 m, (1) 1440 s on the 3.54 m New 
Technology Telescope, (m) 240 s on the 4.2 m William Herschel Telescope, and (n) 600 s on the 2.5 m Isaac 
Newton Telescope. 
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Fig. 2. — Orbital positon plot of active and inactive phases of 176P detailed in Table 1. The Sun is shown 
at the center as a solid dot, with the orbits of Mercury, Venus, Earth, Mars, 176P, and Jupiter (from the 
centre of the plot outwards) are shown as black lines. Solid circles mark positions where 176P was observed 
to be active, while open circles mark positions where 176P was observed to be inactive. Perihelion (P) and 
aphelion (A) positions are also marked with crosses. References: (a) 2005 Oct 24 (Hsieh 2009), (b) 2005 Nov 
26 (Hsieh & Jewitt 2006); (c) 2005 Dec 22-29 (Hsich & Jewitt 2006); (d) 2006 Feb 03-08; (e) 2006 Aug 31 - 
Sep 02; (f) 2006 Dec 11-18; (g) 2007 Jan 27; (h) 2007 Feb 15-16; (i) 2007 Mar 21-22; (j) 2007 May 19; (k) 
2008 Jun 29 - Jul 01; (1) 2009 May 03, where (d)-(l) are from this work. 
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Fig. 3. — Phase functions for 176P. Points are estimated i?-band magnitudes (normalized to heliocentric 
and geocentric distances of 1 AU; tabulated in Table 1) at the mid-point of the full photometric range of the 
nucleus's rotational light curve. Solid circles denote photometry obtained while 176P was visibly active, while 
open circles denote photometry obtained while 176P appeared to be inactive. The dashed line represents a 
least-squares fit (excluding photometry points for which a < 5° where an opposition surge effect is expected) 
to a linear phase function where mfj{l, 1,0) — 15.35 ± 0.08 mag and /3 = 0.038 ± 0.008 mag deg^. The 
solid line represents an lAU {H, G) phase function fit where iJ/j = 15.09 ± 0.05 mag and Gr = 0.15 ± 0.10, 
while the dotted lines indicate the expected range of possible magnitude variations (^ 0.35 mag) due to the 
object's rotation. 



0) 
3 

c 
B 

T3 
C 

I 



O 

m 

< 



14.5 



15 - 



15.5 - 




16 - 



2007 
2007 
o 2007 
• 2007 
A 2007 








0.5 1 

Rotational phase 



1.5 



Fig. 4. — Phase-angle-normalized, reduced-magnitude data {i.e., normalized to a = 0° and R = A = 1 AU) 
for observations of 176P/LINEAR made between 2007 February and May, phased to a rotation period of 
Prot = 22.23 hr. 
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Fig. 5. — Phase-angle-normalized, reduced-magnitude data (i.e., normalized to a = 0° and i? = A = 1 AU) 
for observations of 176P/LINEAR made between 2005 Dec 22 and 2005 Dec 29, phased to a rotation period 



of P., 



22.23 hr. 




Fig. 6. — Simulated lightcurves (solid lines) based on triaxial ellipsoids (shown to the right of each 
lightcurve plot) with axis ratios h/a = 0.54 (red), h/a = 0.56 (green) and h/a = 0.58 (blue) located 
at the orbital configuration of 176P during the 2007 (left panel) and 2005 (right panel) observations. 
The simulations assume obliquity e = 60°, a solstice position vq = 21°. The 2007 and 2005 data 
are overplotted as open and filled circles, respectively. 
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Fig. 7. — Same as Fig. 6 but for solstice position i^q = 0° and ellipsoids with axis ratios b/a = 0.50 
(red), b/a = 0.52 (green), and b/a = 0.54 (blue). 
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Fig. 8. — Same as Fig. 6 but for obliquity e = 45°and an ellipsoid with axis ratio b/a = 0.72 (solid 
orange line). 
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Fig. 9. — Same as Fig. 6 but for obliquity e = 85°, solstice position vq = 40° and ellipsoids with 
axis ratios b/a = 0.33 (red), b/a = 0.40 (green), and b/a = 0.48 (blue). 
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Fig. 10. — Schematic diagram showing how a near-polar jet of directed ejected material can be approximated 
by a wider jet with an effective orientation equivalent to the direction of the rotational pole itself. 




Fig. 11.— Observations of 176P (obtained with Gemini North on (a) 2005 Nov 26 and (b) 2005 Dec 29) 
used to constrain numerical models (§4), where images are shown in the ecliptic coordinate system such that 
the orbital plane of 176P is effectively horizontal and where radiation pressure pushes dust particles to the 
left. 
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Fig. 12. — Images of dust ejection models for 2005 Nov 26 for 10^'^ < /? < 10^^ and different jet directions 
as labeled, where ajet is constant for each row of models and i5jet is constant for each column of models. In 
all panels, the source of emission (i.e., the nucleus) is at the center of each image. 
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Fig. 13. — Images of dust ejection models for 2005 Dec 29 for 10"'^ < /? < 10^^ and different jet directions 
as labeled, where ajet is constant for each row of models and <5jet is constant for each column of models. In 
all panels, the source of emission (i.e., the nucleus) is at the center of each image. 
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Fig. 14. — Comparison of observations (first column of panels) with images of impulsive dust ejection 
models using different particle size distributions (10~^ < f3 < 10^^ in the second column, 10"'^ < P < 10~^ 
in the third column, and 10~^ < /? < 10^^ in the fourth column). Data and models for 2005 November 26, 
2005 December 29, and 2006 February 03 are shown along the first, second, and third rows, respectively. 
All models consist of dust ejected on a single day on November 15 (four weeks after perihelion and 11 days 
before cometary activity was first observed on November 26). In all panels, the source of emission {i.e., the 
nucleus) is at the center of each image. 



